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charges  by  observing  the  detonations  with  a  camera.  The  data  also 
have  interesting  implications  for  consideration  in  theoretical  treat¬ 
ments  of  detonations  in  cylindrical  charges. 


GEORGE  G.  BALL 
Captain,  USN 
Commander 


SOLTR  69-235 


b'" 


v 


4 


TABUS  OP  COm»TS 


IETROD0CTXOI . 

EXPERIMENTAL 

Materials  . 

Pressed  Charge  Preparation  ....  . 

Measurements  . 

Data  Reduction . . . 

RESULTS  AID  DISCUSSIOI 

Detonation  Velocity  as  a  function  of  Density 
Radius  of  Curvature  as  a  function  of  Density 

Variation  of  Rc  with  Charge  Length . 

AP . 

EQ . 

TXT,  p0  -  1.08l  g/cc  . 

Pressed  EDX  . 

Tetxyl  .......  . 

Summary  Rc  vs  £  . 

Estimates  of  Reaction  Zone  Length  . 

SMART  ASD  CQBCWJSIC® . 

REPBRBfCES . 

APPENDIX  . 


Page 

1 


2 

3 

3 

k 

7 

9 

11 

12 

14 

16 

16 

18 

18 

19 

21 

23 


A-l 


ill 


NOLTR  69-235 


TABLE  OF  CONTENTS  (cont’d) 


ILLUSTRATIONS 

Figure  Title  Page 

1  typical  Smear  Camera  Records  of  Curved  Detonation  Fronts  2k 

2  Smoothed  (x,  y)  Data  from  Records  Shown  in  Figure  1  25 

3  Two  Models  for  Expansion  of  Detonation  Front  in  Cylindrical  26 

Charge 

4  D  vs  ^  Curves  for  9  |i  and  25  M  AP  at  d  ~  50  mm  27 

5  Failing  Reactions  of  9  M  AP(N138)  in  Gap  Test  Confinement  28 

6  Rc  vs  fb  Curves  for  9  U  and  25  u  AP  at  d  ~  50  mm  29 

7  Radius  of  Curvature  vs  Length  of  AP  Charges  of  H138  and  30 

N126,  d  =  50.8  mm 

3  Rc  vs  Length  of  Charges  of  NQ-f,  (X5kj)  31 

9  Rc  vs  Length  of  RDX  (x659>  Charges  32 

10  Rc  vs  Length  of  Tetryl  Charges,  d  =  50.8  mm  33 

11  Summary  of  NOL  Radius  of  Curvature  vs  Charge  Length  34 

Measurements 

12  Nominal  Reaction  Zone  Lengths  for  AF  as  a  Function  of  35 

Density  and  Particle  Size 

A1  Extrapolation  of  D  vs  d"1  Data  for  AP(N136),  8.4  u  A- 2 


iv 


KOLER  69-235 


TABLE  OF  CONTESTS  (contfc) 


TABLES 


Table 


Title 


Detonation  Velocity  Data  for  APa  .  J 

Subcritical  Reactions  in  9  ^  AP(N138) . 

OQ 

Radius  of  Curvature  as  a  Function  of  Density  in  AP.  .  3 


Radius  of  Curvature  vs  Charge  Length  for  AP 


Front  Curvature  Found  in  Subdetonation  Reaction  of  , 

AP  Charges . A0 

Radius  of  Curvature  vs  Charge  Length  for  HQ-j!  (X547) 

Radius  of  Curvature  Measurements  on  TWF,  X517  at  . 

„„  >  1.081  g/cc . . te 

Effect  of  Charge  Length  on  Measured  D,  R~  of  RDX,  , 

X659  at  p0  -  1.30  g/cc . ; .  ^3 

Nominal  Reaction  Zone  Lengths  (z)  Computed  from  Eq.  4  ^ 

D  vs  d  Data  for  AP(H136),  8.4  p.  .  A~3 

Aging  Effect  on  AP(N139)>  8.8  u  ..........  A-lf 


NOLTR  69-235 


COMPARISON  OF  CURVATURE  OF  DETONATION  FRONT  IN  AP 
WITH  THAT  POUND  IN  SOME  CONVENTIONAL  EXPLOSIVES 

INTRODUCTION 

Previous  work  has  shown  that  the  radius  of  curvature  (Rq)  of  the  detonation 
front  in  point- initiated  tetryl  (p0  =  1.51  g/cc)  cylinders  increases  with  cylinder 
length  (/)  in  the  manner  to  he  expected  for  spherical  expansion  of  the  front. 

Up  to  the  largest  charge  fired  (H /d  of  4  for  a  =  50.8  mm),  vhere  d  is  the  charge 
diameter,  this  direct  proportionality  between  Rq  and  i  was  observed.  On  the  other 
hand,  Cook  reported2  for  fine-grained,  low  density,  "ideal"  explosives  that  the 
region  of  spherical  expansion  is  followed,  by  a  decaying  rate  of  expansion  until  a 
constant  curvature  is  reached  and  maintained  at  and  above  ~  3- 5-  For 

"non- ideal"  explosives  (e.g. ,  very  coarse,  low  density  TNT  and  coarse  TNT/sodium 
nitrate,  50/50),  the  change  from  spherical  expansion  to  constant  curvature  is 
abrupt  and  occurs  at  i/d  of  1. 5  -  2.0. 

Earlier,  we  developed  a  correction  to  our  measured  detonation  velocities  to 
compensate  for  the  curvature  of  the  initiating  shock  front. 3  This  correction  was 
derived  on  the  assumption  of  uninterrupted  spherical  expansion  of  the  detonation 
front.  Consequently  the  validity  of  this  correction  depends  upon  its  application 
only  to  explosives  exhibiting  such  behavior.  Thus,  one  objective  of  the  present 
study  was  to  examine  Rc  vs  f/d  data  of  several  explosives  to  determine  which  ones 
exhibit  the  spherically  expanding  fronts. 

Another  objective  was  to  obtain  the  detailed  information  about  curvature  of 
the  detonation  front  needed  to  inteipret  correctly  (as  particle  velocity)  the 
measurements  obtained  by  the  electromagnetic  method.  This  method  is  in  process 
of  development;  it  depends  on  a  thin  metal  foil  moving  at  the  material  velocity 
of  the  front.  Best  results  (no  correction)  will  be  obtained  for  planar  detonation 
fronts.  Acceptable  results  can  be  obtained  by  correcting  the  measurements  accord¬ 
ing  to  independently  measured  Rq.  The  curvature  must  be  known,  must  be  reproducible, 
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and  must  be  small  enough  to  set  a  practical  length  of  foil  (e.g.,  5  -  10  mm)  into 
uniform  motion  within  a  very  small  delay  period  (e.g.,  30  nsec. ).  In  particular 
we  wished  to  know  Hq  of  ammonium  perchlorate  (AP)  as  a  function  of  particle  size 
(5),  loading  density  (p6),  and  charge  length  so  that  the  electromagnetic  method 
can  be  used  to  study  what  is  generally  considered  to  be  the  relatively  long  re¬ 
action  zone  of  this  explosive. 

In  addition  to  the  primary  objectives  of  obtaining  data  necessary  for  the 
electromagnetic  method  and  for  valid  correction.,  of  our  detonation  measurements 
we  irished  to  examine  any  appropriate  curvature  data  in  the  theories*1^  Which  relate 
Rc  to  reaction  zone  length,  a.  Pew  applications  of  these  theories  have  been  made 
because  appropriate  Rc  data  are  not  usually  available. 

EXPERIMENTAL 

Materials 

Five  batches  of  AP  were  used  in  this  work}  all  were  of  the  propellant  grade 
and  contained  0.2  -  0. 5$  tricalcium  phopphate  (TCP).  Four  batches,  freshly  ground 
on  the  same  schedule  at  the  Naval  Ordnance  Station,  Indian  Head,  Maryland,  had  a 
nominal  average  particle  size  of  9n  by  micromerograph.  These  were  N136  (8.4m), 

N138  (8.  9  m),  N139  (8*3  m)  and  Nl4o  (8.3  m)«  The  fifth  batch,  N126  was  originally 
assigned  an  average  particle  size  of  25  M°  when  it  was  acquired  several  years  ago. 
After  putting  a  sample  of  this  material  through  a  coarse  screen  to  remove  a  few 
large  agglomerates,  NOS  reports  that  its  micromerograph  value  is  still  25  m • 

Three  organic  H.E.  were  used  in  this  work  and  all  were  procured  under  the 
respective  military  specifications.  RDX,  X659>  was  ordered  as  Type  B  Class  E 
(98$  through  a  No.  325  screen)  Holston  Production;  an  analogous  lot  of  HMX  had 
a  weight  mean  particle  size  of  14  m«  TNT,  X517>  is  a  coarse  H.E.  with  a  mean 
5  of  200  m;  see  Table  7.  Nitroguanidine  (NQ)  X5^7  is  the  low  bulk  density  form 
described  and  used  in  earlier  work. ^ 
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Pressed  Charge  Preparation 

The  explosives  were  handled  and  made  into  cylindrical  charges  as  in  the  pre¬ 
vious  work®  except  that  AP  charges  were  given  special  treatment  to  avoid  moisture 
pick-up.  In  preparing  these  charges,  AP  was  taken  directly  from  moisture- 
vapor  proof  storage  bags  and  loaded  directly  into  the  mold;  if  for  any  reason 
AP  required  drying,  it  was  dried  in  a  vacuum  oven.  Any  AP  charges  not  fired  as 
soon  as  possible  after  preparation  were  packaged  in  Saran  and  stored  at  30°C  until 
they  could  be  fired. 

In  the  present  work,  the  charge  diameter  was  kept  at  50  -  51  nun,  but  both 
length  and  the  loading  density  were  varied.  The  total  charge  train  was  made  up 
of  a  detonator  (Engineers  Special  Hercules  J-2),  a  50/50,  1-56  g/cc  pentollte 
booster  (diam.  of  main  charge,  25.4  mm  long),  and  the  main  charge.  In  detonation 
velocity  measurements,  the  main  charge  was  generally  followed  by  a  second  booster 
pellet  to  serve  as  an  explosive  witness. 

Measurements 

A  70  ram  smear  camera  (writing  speed  about  4  mm/// sec)  was  used  to  measure 

both  the  detonation  velocity  (d)  and  the  radius  of  curvature.  All  records  were 

made  with  Tri-X  film  developed  in  Acufine. 

For  measurements  of  D,  the  charge  was  mounted  vertically  with  the  detonator 

at  the  top.  In  a  veiy  few  cases,  Magic  Tape  was  used  as  a  flasher*.  The  usual 

8 

procedure  and  typical  records  are  described  in  previous  work. 

For  measurements  of  Rq,  the  charge  was  mounted  horizontally  on  a  Styrofoam 
V-block.  The  camera  slit  was  focused  on  the  end  face  of  the  main  charge,  and  was 
perpendicular  to  the  charge  axis.  The  maximum  writing  speed  of  4  mm///  sec  was 

*  If  possible,  a  flasher  is  avoided.  For  granular  charges,  especially  those 
with  chalky  surfaces,  it  is  difficult  to  produce  an  air  gap  of  uniform  thick¬ 
ness  and  hence  a  record  of  uniform  luminosity. 
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used  but,  even  so,  the  horizontal  (time)  c oniponent  of  the  trace  varied  only  from 
2  to  11  mm  on  the  film.  Obviously  a  higher  writing  speed  would  improve  this 
measurement.  Good  traces  were  obtained  Without  using  a  flasher.  Pig.  1  shows 
two  typical  records. 

Data  Reduction 

Reading  and  reduction  of  the  smear  camera  velocity  records  is  briefly  des¬ 
cribed  in  Ref.  8.  Details  of  using  the  Universal  Telereader  and  auxiliary  equip¬ 
ment  are  set  ourt  in  Ref.  9*  Thirty  to  40  points  were  read  on  each  record  for  the 
work  reported  here.  For  the  most  part,  the  distance-time  data  were  fitted  to  a 
straight  line  (by  a  least  squares  procedure)  in  order  to  determine  the  detonation 
velocity.  The  use  of  the  straight  line  was  justified  by  the  adequacy  of  the  fit. 

For  short  charges,  a  quadratic  was  fitted  to  the  data,  also  by  least  squares.  The 
velocity  at  any  given  point  on  the  observes,  part  of  the  charge  wa3  then  determined 
by  evaluating  the  derivative  of  the  quadratic. 

The  wave  front  records  were  read  on  the  same  machine  as  the  velocity  records. 
The  cross  wires  are  set  on  the  streak  camera  trace  as  shown  in  Pig.  lb.  That  is, 
the  vertical  wire  is  set  on  the  left-most  part  of  the  trace,  so  it  coincides  with 
line  sib .  Then  the  horizontal  wire  is  moved  up  to  a  point  where  the  trace  is 
distinct,  point  p.  Thus  point  p  is  the  origin  for  the  readings,  and  about  ’+0  points 
are  read  across  the  trace.  Any  tilt  of  the  slit  with  respect  to  the  axis  of 
rotation  of  the  camera  mirror  is  compensated  for  by  reading  the  coordinates  of 
the  ends  of  the  image  of  the  slit,  ef  and  gh. 

The  ordinates  of  the  upper  and  lower  edges  of  the  image  of  the  end  of  the 
charge  are  also  recorded.  The  horizontal  wire  is  set  successively  on  ij  and  ET. 
These  latter  are  used  to  locate  the  center  of  the  charge  so  that  the  origin  for 
the  trace  can  be  moved  to  point  o,  see  pig.  lb,  and  the  line  ot  is  the  time  axis. 

The  scale  in  the  picture  is  read  to  give  the  proportionality  between  distance  on 
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the  charge  to  distance  on  the  screen  of  the  Telereader.  The  writing  speed  of  the 
camera  and  the  magnification  of  the  film  in  the  Telereader  give  the  factor  required 
to  compute  time  from  the  values  of  the  abscissa  as  read. 

In  previous  work  at  N011,  data  similar  to  those  described  above  were  found 

to  represent  a  spherically  expanding  front;  in  cross  section,  the  fronts  were 

2 

circular.  Also,  Cook  found  that  short  charges  of  numerous  explosives 

gave  spherically  expanding  fronts.  For  longer  charges,  however,  Cook  found  that 
although  the  fronts  were  spherical,  the  radius  of  curvature  did  not  increase  with 
charge  length.  Thus  it  was  expected  that  the  data  from  the  present  work  could  be 
represented  by  a  circle  (because  we  observe  a  cross  section).  This  was  confirmed 
by  the  adequacy  of  the  corresponding  fit  obtained.  See  Fig.  2  for  the  results 
for  the  two  shots  represented  by  the  photographs  in  Fig.  la  and  lb. 

When  the  front  shows  continuous  spherical  expansion,  the  curvature  is  un¬ 
steady  in  the  sense  that  it  changes  with  additional  charge  length.  The  change 
with  length  will,  however,  decrease  with  increasing  charge  le*  th.  In  contrast 
to  this  behavior,  some  charges  show  a  spherical  expansion  of  the  front  only  for 
short  distances,  after  which  Rc  becomes  constant  and  independent  of  length  of 
travel.  It  is  important  to  be  able  to  recognize  these  two  different  modes  of 
propagation.  This  can  be  done  by  superimposing  streak  camera  records  which  have 
been  taken  on  charges  of  different  length  (after  minor  corrections  for  tilt  and 
differences  in  magnification  have  been  made).  If  the  records  coincide,  the  front 
is  steady  for  that  material  and  geometry. 

Another  way  of  determining  if  the  shape  of  the  detonation  front  is  constant 
is  to  observe  the  velocity  of  detonation  with  the  streak  camera.  It  has  been 
shown''"0  that  for  a  spherically  expanding  front,  this  method  gives  an  apparent 
velocity  which  is  greater  than  the  actual  detonation  velocity.  An  alternate 
statement  is  that  when  the  velocity  along  the  periphery  of  a  cylinder  is  not 
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constant,  the  front  may  be  expanding  spherically.  Thus  careful  measurements  of 
the  velocity  along  the  periphery  of  a  charge  can,  in  principle,  be  used  to  dis¬ 
tinguish  between  steady  and  unsteady  curvature  of  the  detonation  front.  We  may 
not  have  sufficient  resolution  to  see  any  change  in  D  with  spherical  expansion; 
in  most  cases,  D  appeared  to  be  constant  within  our  limits  of  resolution. 

The  radius  of  curvature  is  computed  in  different  ways  for  the  two  cases 
discussed  above.  The  radial  expansion  case  is  diagrammed  in  Fig.  3a  where  the 
detonation  front  has  just  reached  the  end  of  the  charge.  The  radius  of  curvature 
is  Rc,  the  distances  sF  and  ac.  It  is  assumed  that  the  detonation  velocity,  D, 
is  known,  and  constant,  so  that  cd  is  D  x  t,  where  t  is  the  time  measured  from 
the  appearance  of  the  front  at  point  b.  This  time  at  off-axis  point  d,tg,  and 
the  abscissa,  xd  are  known  from  the  streak  camera  record.  Thus 

Rc2  +  xd2  a  (Rc  +  D  x  *d)2  (1) 

so  that  Rc  could  possibly  be  computed  from  only  one  point  on  the  record. 

More  representative  results  are  obtained  by  using  the  method  of  least  squares 
so  that  the  subscript  d  takes  on  all  integer  values  from  1  to  the  total  number 
of  points  read  on  the  film. 

Note  that  in  Fig.  3»,  the  front  propagates  in  the  direction  of  its  normal. 
For  the  steady  case,  see  Fig.  3b,  the  front  moves  as  a  rigid  body  so  that  cd  is 
D  x  t  and  the  path  is  parallel  to  the  axis  of  the  charge.  Fof  this  case  the  co¬ 
ordinates  of  a  point  on  the  front,  such  as  c,  are  given  by  (x&,-  D  x  td),  where 
both  x,j  and  t^  represent  a  particular  point  on  the  trace.  All  data  points  are 
treated  in  this  way  so  that  a  new  set  of  data  is  generated  which  represents  the 
arc  shown  in  Fig.  3h.  These  are  used  to  evaluate  h,  k,  and  Rc  in  the  equation 
of  a  circle, 

(x  -  h)2  +  (y  -  k)2  =  Rc  (2) 
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by  least  squares.  T'-’is  generalized  equation  per  i:‘.ts  us  to  obtain  better  fits  to 
data  from  charges  w  *  c,  for  one  reason  or  another,  are  not  perfectly  uniform  or 
perfectly  aligned  in  the  optical  system.  It  also  helps  to  compensate  for  small 
asymmetries  which,  can  be  introduced  by  the  detonators,  i.e.,  these  only  ay...TXi- 
mate  point  initiation.  Pig.  2  shows  the  results  obtained  from  the  records  s’. own 
in  Pig.  1  for  AP  charges  which  exhibit  a  constant  Rc.  In  this  figure,  t..e  solid 
line  trace  is  for  data  smoothed  to  fit  Eq.  (2).  The  points  shown  are  derived 
(see  above)  from  the  actual  readings  of  the  records. 

The  adequacy  of  the  fit  for  either  Eq.  (l)  or  (2)  is  determined  by  examining 
plots  such  as  Pigs.  2a  and  2b,  and  the  computed  q.m. e.  (quadratic  mean  error), 
which  is 


i.m.e.  = 


n  .  ,2 

2  "  >'ic) 


AK 

J  ./<n 


P) 


(3) 


where  y^  =  observed  value  of  ordinate 

yic=  value  of  ordinate  computed  from  the  fitted  equatiot. 
n  •=  number  of  data  points 

p  =  number  of  parameters  in  the  equation  to  be  fitted 
Values  of  the  q.m. e.  have  all  been  in  an  acceptable  range;  for  example,  the 
values  for  the  fits  shown  in  Fig.  2  are  0.13  mm  for  values  of  y  which  vary  from 
-25  to  +25  mm.  The  computer  also  gives  an  estimate  of  the  standard  deviation  of 
each  of  the  parameters  being  estimated1^.  These  have  been  converted  to  percentages 
and  are  reported  in  the  tables  listing  the  results. 

RESULTS  AND  DISCUSSION 


Detonation  Velocity  as  a  Rmction  of  Density 

D  vs  pq  curves  were  obtained  on  the  9  and  25 fi  AP  at  d  ~  50  mm,  the  constant 
diameter  selected  for  studying  the  curvature  of  the  front.  Data  for  the  9m  mater¬ 
ial  were  obtained  on  batches  NI36  and  N138,  both  ground  on  the  same  schedule  at 
NOS;  these  two  batches  are  assumed  almost  indistinguishable  in  behavior.  Data 
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for  the  25m  AP  were  obtained  on  U126  for  which  a  curve  had  been  obtained  several 
years  ago.^  In  view  of  some  evidence  that  the  behavior  of  AP  changes  with  age3^, 
it  was  felt  that  a  second  determination  was  necessary.  The  data  obtained  are 
given  in  Table  1  and  plotted  in  Jig.  4.  These  D  values  have  not  been  corrected 
for  non-planar  initiation  because,  as  the  next  section  will  show,  such  a  correc¬ 
tion  is  not  valid  for  these  AP  charges. 

In  similar  studies  of  other  H.E. ,  in  particular  the  method  of  charge 

preparation  has  seemed  to  affect  some  of  the  results.  Consequently,  the  method 
of  conduction  (by  hand,  by  hydraulic  press,  or  by  isostatic  press)  has  been  indi¬ 
cated  on  the  graphs  as  well  as  in  the  tables.  The  trends  (D  vs  pq )  of  Fig.  4 
seem  continuous  despite  changes  from  one  method  of  charge  preparation  to  another. 
The  curve  for  the  fine  AP  parallels  the  ideal  curve  in  the  range  0. 6  to  1 . 2  g/ cc 
after  which  it  shews  a  moderate  bending  until  it  ends.  The  charge  at  1.4l  g/cc 
failed,  but  this  does  not  necessarily  mean  that  the  critical  density  has  been 
exceeded.  In  this  work  we  are  using  only  half  the  booster  length  used  in  earlier 

O 

work.  Hence  the  failure  could  result  from  inadequate  boostering  or  from  exceed¬ 
ing  the  critical  density.  The  farmer  seems  more  probable  because  the  present 
9m  AP  D  vs  p0  curve  (d  ~  50  mm)  almost  coincides  with  the  10m  AP  curve  at 
d  *  76  mm,  pQ  <1.3  g/cc  in  the  previous  work.  This  is  an  indication  that  the 
9m  AP  is  finer  than  the  10m  AP  as,  of  course,  it  should  be.  For  the  same  reason 
its  detonabllity  curve  should  lie  at  lower  values  of  nc  and  higher  values  of  #c 
than  the  curve  of  the  less  fine  material.  Finally,  the  10m  AP  at  d  *  50.8  mm, 
had  ~  l.V?5  g/cc;  consequently,  the  9m  AP  should  have  Pc  >1.48  g/cc  at 
d  —  50  mm. 

It  is  not  important  for  the  present  work  to  determine  P c  for  this  9m  AP  at 
d  ~  50  ram,  and  no  effort  has  been  made  to  do  so.  However,  it  is  of  interest  that 


8 


HOME  69-235 


in  the  gap  test  confinement,  this  AP  detonates  at  X. 584  g/cc^  and  fails  (is  dead 
pressed)  at  1.590  g/cc.  The  critical  density  in  the  confinement  of  the  gap  test 
must,  therefore,  be  quite  near  1.59  g/cc.  These  limits  were  determined  with  pin 
measurements  which  resulted  in  two  good  examples  of  the  failing  reaction  obtained 
in  strongly  boostered  AP  when  experimental  conditions  make  it  sub-critical. 

Table  2  and  Fig.  5  ehcv  the  fading  reactions  observed  in  AP  (H138). 

The  curve  for  the  25m  AP  is  very  similar  to  that  obtained  when  this  mater¬ 
ial  was  first  studied.  In  fact,  it  is  identical  at  and  to  the  right  (toward 
higher  p0 )  of  the  maximum  in  D.  The  difference  appears  at  lower  densities  and 
amounts,  at  most,  to  0. l4  mm/uaec.  Several  years  of  aging  has  reduced  D  at  the 
lowest  pQ  (0.9  g/cc)  by  this  amount.  These  data  and  one  other  example  of  an 
aging  effect  In  a  finer  AP,  given  in  the  Appendix,  confirm  the  aging  effect  sug¬ 
gested  in  Appendix  B  of  Ref.  (12).  The  particle  size  effect  on  D  is  the  same 
order  of  magnitude  as  that  established  in  previous  work. 0 
Radius  of  Curvature  as  a  Function  of  Density 

Table  3  contains  the  data  Rg  vs  p0  obtained  for  the  9  fi  and  25  p  AP.  The 
radius  of  curvature  was  computed  in  two  ways;  (l)  fitting  all  the  x-y  data 
across  the  charge  surface  (d  —  50  mm)  to  Eq.  2,  and  (2)  fitting  only  half  the 
data  covering  the  inner  portion  of  the  charge  surface  (d  ~  25  ram)  to  Eq.  2.  This 
was  done  because  the  quite  commonly  used  curved  front  theory  of  Eyring  ot  al^, 
proposes  a  large  edge  effect.  In  fact,  the  effect  in  the  theory  as  originally 
formulated  is  so  large  that  its  presence  would  make  a  fit  of  all  data  to  either 
equation  meaningless. 

Exam"  jation  of  the  data  of  Table  3  shows  that  it  is  quite  possible  to  fit 
most  of  the  x-y  data  over  the  entire  surface  with  a  percent  error  of  1.6$  Rc  or 
less.  The  analogous  fit  for  the  central  part  of  the  same  surface  does  give  a 
soiuewhat  greater  Rc  which  way  indicate  a  small  edge  effect.  (Many  records, 
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typified  by  Pig.  2,  show  a  small  but  definite  increase  in  the  curvature  of  the 
front  as  it  approaches  the  charge  surface).  The  percent  error,  in  this  case,  is 
somewhat  larger  than  that  for  treating  the  whole  surface,  but  generally  only  to 
the  extent  that  would  be  expected  by  halving  the  number  of  measurements  treated. 

If  the  large  edge  effect  exists,  it  cannot  be  resolved  on  our  records.  Rg  deter¬ 
mined  by  either  method  showed  the  same  trends.  Hence  we  chose  to  work  here  with 
that  determined  across  the  entire  charge  surface  with  the  advantage  of  its  smaller 
percent  error.  Only  values  so  determined  are  tabulated  and  plotted  in  the  rest  of 
the  report. 

The  change  in  curvature  noted  above  is  not  consistent  with  the  Wood  &  Kirkwood^ 
theory  nor  with  the  Ryring5  theory.  However,  these  data  may  be  useful  in  modifying 
the  theories  as  time  permits.  For  example,  a  condition  can  be  put  on  the  way  the 
detonation  front  meets  the  edge  of  the  charge  in  the  Fyring  theory.  This  results 
in  a  major  reduction  of  the  confuted  curvature.  It  is  not  known  if  the  theory  of 
Wood  and  Kirkwood  can  be  modified  to  give  radii  of  curvature  whose  values  differ 
at  the  edge  of  the  charge  from  that  at  the  center  of  the  charge.  For  these  reasons, 
the  radii  of  curvature  over  the  central  half  of  the  AP  charge  faces,  given  in 
Table  Sjinay  be  useful  in  testing  any  future  modification  of  the  curvature  theories. 

Fig.  6  shows  the  plot  of  R,,  vs  p0  for  the  9  n  and  25  /u  AP.  The  small  varia¬ 
tions  in  diameter  (2$  or  less)  have  been  ignored  on  the  assunption  that  any  corres¬ 
ponding  effect  on  Rc  would  be  undetectable.  No  distinction  has  been  made  in 
charge  lengths  of  203  to  457  mm  because,  as  data  in  the  next  section  will  she**, 
charges  from  both  batches  of  AP  attain  and  maintain  a  constant  Rq  over  this  range 
of  length.  As  in  the  case  of  Fig.  4  (D  vs  pQ),  the  method  of  charge  preparation 
has  been  indicated  in  Fig.  6  It  seems  to  have  more  effect  on  Rc  than  on  D. 

Over  the  range  of  loading  density  for  which  charges  were  compacted  in  the 
isostatic  press,  the  trend  shown  in  Fig.  6  for  both  batches  of  AP  is  a  decrease 
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in  Rc  (increase  in  curvature)  with  an  increase  in  Pq.  This  is  opposite  in  direc¬ 
tion  to  the  trend  Cook^  reported  for  "ideal"  explosives;  he  gives  Rg  values  for 
the  "non-ideal"  at  only  one  density.  Fig.  6  shows  a  particle  size  effect  on  Rc; 
for  charges  ccopacted  in  the  isostatic  press,  Rc  decreases  with  increase  in 
initial  particle  size.  This  is  the  trend  reported  by  Cook  for  both  ideal  and 
non- ideal  explosives. 

For  the  finer  AP  (N133)*  the  trend  found  with  the  charges  prepared  in  the 
isostatic  press  can  be  reasonably  extrapolated  to  the  p0  range  of  the  charges 
prepared  in  the  hydraulic  press  (See  dashed  line  of  Fig.  6).  The  scatter  of  the 
data  in  the  latter  region  seems  a  bit  more  than  twice  that  of  the  former.  This 
agrees  with  the  poorer  reproducibility  found  in  hydraulically  pressed  charges. 

For  the  coarser  AP  (IT126),  no  analogously  reasonable  extrapolation  is  available, 
on  going  from  charges  prepared  in  the  isostatic  press  to  those  which  were  hand- 
packed.  For  this  material,  loose  packing  seems  to  overwhelm  the  density  effect 
shown  by  pressed  charges. 

Within  the  higher  density  region  (1.19  to  1.38  g/cc)  where  the  data  seem  to 

establish  valid  trends,  there  is  no  apparent  con elation  between  R,  and  D.  In 

particular,  Rg  does  not,  like  D,  show  a  change  indicating  an  approach  to  failure 

conditions  nor  does  it  extrapolate  to  a  value  as  low  as  0.5  d  (25  mm)  at  pc  the  Rq 
o 

reported  at  dc. 

Variation  of  Re  with  Charge  Length 

Berbaps  the  most  interesting  result  of  the  present  study  is  the  variation 
of  Rc  with  charge  length  found  for  several  explosives.  These  will  be  described 
in  turn. 
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AP 

Table  4  contains  radii  of  curvature  for  both  9  u  and  25  u  AP  in  charges 
varying  in  length,  /,  from  2.5  to  23.0  cm.  Values  of  Rc  in  the  fifth  column 
are  determined  from  fitting  Eq.  2  to  points  all  the  iray  across  the  charge. 

For  H  <12. T  cm,  Eq.  1  has  also  been  used  in  the  same  manner,  the  results  being 
reported  in  column  6.  For  these  data,  there  is  little  difference  between  the 
results  obtained  from  the  two  equations.  Because  the  AP  charges  very  rapidly 
develop  a  constant  l^,  in  the  following  discussions  we  will  only  refer  to  the 
results  obtained  from  Eq.  2  (col.  5)j  these  are  plotted  in  Fig.  7.  Discussion 
of  this  figure  will  be  more  understandable  if  we  first  review  briefly  some  of 
our  general  knowledge  of  shock- to-detonation  transition. 

The  shock  initiation  of  detonation  depends  upon  the  area  of  shock  front  in 
the  acceptor  as  well  as  the  pressure- time  history  which  that  shocked  area  experi- 
ences.  According  to  the  particular  combination  of  these  factors,  the  detonation 
induced  in  the  acceptor  nay  be  temporarily  under  or  overdriven.  In  the  first  case, 
the  shock  induced  reaction  is  highly  exothermic  but  is  sub-detonation.  Consequently, 
the  energy  of  reaction  exceeds  the  energy  lost  by  any  dissipative  processes  until 
the  reaction  has  built  up  to  its  steady  state.  In  the  case  of  overdrive,  a  strong 
shock  initiates  reaction  at  a  higher  temperature  and  pressure  than  the  steady  state 
C-J  values.  Consequently,  the  excess  energy  is  dissipated  until  the  steady  state 
condition  is  reached.  Overdrive  (or  underdrive)  decreases  to  a  negligible  amount 
by  the  time  the  disturbance  has  travelled  to  (i/d)  =  2  even  when  the  donor  and 
acceptor  are  badly  mismatched.3^ 

The  data  of  Table  4  show  that,  for  AP,  Rc  decreases  to  the  constant  value 
each  attains  at  (£/ d)  >  1«5«  A  similar  decrease  in  D  can  be  detected  at  the  lowest 
density.  Thus  our  initiation  system  of  a  detonator  followed  by  a  25.4  mm  long 
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pentolite  booster- overdrives  the  AP  detonation,  and  the  overdrive  increases  with 
increasing  particle  size.  The  detonation  front  in  AP  does  not  exhibit  a  spherical 

p 

expansion,  but  rapidly  attains  constant  a  behavior  reported  by  Cookr  for  non¬ 
ideal  explosives.  At  small  values  of  i,  however,  the  Ref.  (2)  curves  followed 
the  spherical  expansion  line,  i.e.,  detonator  initiation*  (approximately,  point 
initiation)  resulted  in  an  underdriven  detonation  which  built  up  to  the  steady 
state# 

It  is  concluded  from  the  constant  plateau  values  found  in  each  case  that  plane 
rave  propagation  of  steady-state  detonation  in  AP  is  impossible  in  50  mm  diam. 
charges.  Thus  a  plane  wave  booster  would  be  expected  to  affect  the  wave  shape  only 
in  the  transitional  region  in  which  the  over  boostering  is  fading  out.  This  was 
confirmed  by  shots  772  and  773  in  which  ca.  9  M  AP  (RlUo)  at  P0  =  1.29  g/cc  was 
initiated  with  a  PWB.  At  (i/d)  of  1.0  and  2.0,  the  Rc  values  were  6l  and  66  mm, 
respectively.  The  latter  value  falls  on  Mg.  7  about  where  it  should  for  9  M  AP, 
i.e.,  linear  interpolation  on  the  density  to  1.29  g/cc  gives  an  Rc  value  for  AP 
(N138 )  of  65  mm.  Thus  the  plateau  region  of  Mg.  7  seems  independent  of  the 
original  source  of  initiation.  Hence  there  is  no  point  to  initiating  AP  with  a 
PWB  for  the  electromagnetic  studies.  On  the  other  hand,  it  seems  safe  to  assume 
that  constant  D  and  Fcc  will  have  been  attained  at  (//d )  *1.5  instead  of  only  at 
(i/a)  a  3.0  as  we  have  been  assuming.  For  steady  state  values  then,  we  can  use 
shorter  charges  than  we  had  thought  necessary. 

In  the  corrections  we  have  recently  used  on  our  detonation  velocity  measure¬ 
ments,  we  assumed:  (l)  spherical  expansion  of  the  detonation  front  and  (2)  an 

*  Where  a.  booster  was  necessary,  the  detonator  vas  placed  in  a  19  mm  deep  well  in 
a  25A  ran  booster.^  Unfortunately,  Cook  did  not  indicate  whether  his  coarse  THT 
and  THT/k'aNOj  mix  were  cap  sensitive  or  not.  If  not,  the  detonator  would  still 
have  been  only  6.4  mm  from  the  main  charge  surface  and  hence  approximated  point 
initiation  better  than  the  initiator  of  the  present  work. 
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increase  in  as  detonation  passed  from  the  "booster  (donor)  to  the  main  AP  charge 
Receptor).  Fig  7  shows  that  both  of  these  assumptions  are  wrong  for  AP.  For  this 
explosive,  the  uncorrected  velocity  (measured  optically  or  with  pins  along  the 
charge  surface)  is  also  the  axial  velocity. 

Finally,  several  non-detonating  charges  (fading  r« action)  were  examined. 

These  data  in  Table  5  show  no  obvious  change  in  Rc  witt  failure  to  detonate; 
the  x-y  data  can  be  fitted  to  Eq.  1  or  Eq.  2  as  well  as  those  from  detonating 
charges.  However,  the  luminosity  of  the  trace  shows  interesting  changes  with 

charge  length.  At  /*  25.4  mm  or  (//d)  ®  0. 5,  the  effect  of  the  booster  is  domi¬ 

nant  here  as  it  is  for  the  detonating  charges.  Rq  is  71  ram  (Shot  75*0  as  compared 
to  72-7**  mm  in  Fig.  7-  In  other  words,  at  (//d)  *  0.5,  Rc  is  the  same  regardless 
of  particle  size,  loading  density,  and  whether  the  disturbance  is  a  strong  shock 
or  a  true  detonation.  At  (/,/d)  =  1.5  for  the  26  n  AP  i.t  «  l.**35>  Rc  is  still 

in  the  range  of  the  Fig.  7  curves  and  the  wave  profile  in  AP  at  this  density  and 

particle  size  has  faded  out  at  both  edges  (See  Fig.  2,  Shot  731)  so  that  only  0.9 
d  is  covered.  An  even  better  example  of  the  progress  of  fading  of  the  reaction 
from  the  surface  of  the  charge  to  its  central  axis  is  shown  by  this  AP  at 
P0  =  1.40  -  1.4l  g/cc.  This  density  is  nearer  the  failure  point  (ca.  1.385  g/cc, 
see  Fig.  4)  and  consequently  the  failing  reaction  persists  longer.  At  (// d)  =  3, 
the  wave  profile  still  covers  the  diameter,  but  by  (£/ d)  =  3.5  only  a  very  faint 
trace  can  be  seen  which  covers  just  the  central  third  of  the  diameter. 

NQ 

The  second  explosive  studied  was  low  bulk  density  nitroguanidine,  NQ-/,  X5**7. 
This  was  studied  at  56$  and  92$  OMD  rather  than  at  60  and  69$  USD  as  was  the  case 
for  AP.^  This  greater  spread  in  density  may  make  the  results  differ  significantly 
from  those  for  the  AP.  Data  on  results  for  both  densities  of  NQ-Z  are  given  in 
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Table  6;  the  radii  of  curvature  are  shown  in  Pig.  8.  Eq.  1  was  used  to  fit  the 
data  for  the  high  density  charges,  Eq.  2  for  the  low  density  charges*.  In  other 
words,  the  shock  from  the  1. 65  g/ cc  material  seems  to  expand  more  nearly  radially 
than  does  the  1. 0  g/ cc  material.  The  solid  diagonal  line  in  Pig.  8  represents 
simple  radial  expansion  with  no  regard  for  the  booster  charge.  Hence  it  is  not 
surprising  that  the  point  for  j?  =  76.2  mm  lies  above  the  solid  line.  The  dashed 
curve  passes  through  two  points  as  a  straight  line  and  extrapolates  to  an  effec¬ 
tive  initiator  length  of  22  mm  (i.e.,  1=  -22.0  for  Rc  =  0),  which  is  close  to 
the  25. H  mm  booster  length.  This  is  reasonable  since  the  pentolite  booster  and 
the  higher  density  HQ  have  nearly  the  same  impedance.  Radial  expansion  apparently 
breaks  down  somewhere  between  l  =  127  and  229  mm,  or  between  (i/d)  *  2.5  and  4. 5. 
If  the  actual  booster  length  is  added  to  the  length  of  the  charge,  we  have  break¬ 
down  between  3*0  and  5.0  for  /*■  i+  25*  ^  mm.  These  numbers  are  more  sug¬ 

gestive  than  accurate.  For  example,  note  the  error  bars  which  have  been  placed 
on  the  points  for  pQ  =  1.65  g/cc.  These  errors  are  based  on  one  standard  devia¬ 
tion  of  the  computed  value  of  the  quantity  Rq  in  Eq.  1.  Obviously  we  need  more 
data  before  we  can  be  certain  about  the  behavior  of  this  material. 

The  behavior  of  NQ ,-£  at  a  density  of  1. 0  g/cc  is  much  more  similar  to  that 
of  AP.  Equation  2  was  used  to  fit  these  data  and  the  results  are  shown  in  Fig. 8 
also.  We  find  a  constant  value  of  Rc  of  mm  beyond  ji  =  76. 2  mm.  Hence  this 
material  has  a  constant-  Rc  for  (£/d)  >1.5*  As  in  the  case  of  the  AP,  this  mater¬ 
ial  is  possibly  overboostered.  Since  our  data  show  only  the  constant  value  of  Rc, 
the  estimate  of  effective  initiator  length  (l8  mm)  can  be  only  a  lower  limit 

*  Data  'of  Table  6  show  that  fitting  the  data  for  1.65  g/cc  with  Eq.  2  give? a 
trend  of  Rq  vs  £/&  that  cannot  be  readily  interpreted.  The  use  of  Eq.  1,  however, 
leads  to  the  interpretation  of  the  text. 
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approximation.  Again,  more  data  are  needed  to  clarify  the  situation  completely. 

We  have,  however,  enough  data  to  guide  us  in  future  work  on  1.0  g/cc  density 
charges  of  HQ-/  in  the  electromagnetic  experiments. 

It  is  evident  that  for  HQ-/,  Rq  increases  as  pQ  increases,  the  same  trend  as 
that  reported  for  ideal  explosives  by  Cock^  and  the  reverse  of  the  trend  found  for 
AP.  The  trend  with  density  indicates  that  HQ-/  at  pQ  *0.4  g/cc  will  have 
Rc  <  94  ram  ?»nd  constant  Rc  at  (//d)  >1.5.  Hence  the  successful  use  of  such  low 
density  HQ  in  a  plane  wave  generator1^  must  have  resulted  from  overboostering 
the  HQ. 

y  »  1-Oft!-  g/cc. 

Because  HQ  seems,  in  seme  ways,  an  atypical  organic  H.E. ,  radius  of  curvature 
measurements  were  also  carried  out  on  low  density  TKT,  X517.  The  data  are  given 
in  Table  7.  Although  the  scatter  in  Rc  is  large,  this  would  be  expected  because 
of  three  difficulties  in  controlling  charge  preparation.  These  are  large  particle 
size  (ca.  200  p. ),  low  charge  density  (-^  1.0  g/cc),  and  compaction  in  the  hydraulic 
press,  all  factors  which  favor  local  irregularities  in  the  charge.  With  allowance 
for  the  larger  scatter,  it  seems  clear  that  low  density  TOT?  behaves  like  AP,  i.e., 
has  a  constant  Eg  at  (£/d)  >  1.5. 

Pressed  KDX 

The  lowest  density  charges  of  KDX  which  could  be  prepared  in  the  isostatic 
press  were  those  at  about  1.3  g/cc.  This  material, X659>  is  fine  grained  and  the 
charges  were  expected  to  be  of  good  quality.  That  is,  isostatic  pressing  of  a  fine, 
chemically  homogeneous  material  usually  results  in  uniform  charges  as  was  true  in 
this  case.  Data  for  the  experiments  are  given  in  Table  8;  these  are  for  detona¬ 
tion  velocity  in  KDX  at  pQ  =1.3  g/cc  and  for  curvature  studies  in  both  boostered 
charges  of  KDX  and  detonator  initiated  charges.  Equation  1  was  applied  to  data 
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from  these  shots  and  the  resulting  values  of  Rq  are  shown  in  Fig.  9.  Again  the 
error  bars  are  given  by  the  computer  program  which  fits  the  data  —  we  have  not 
analyzed  the  situation  for  experimental  errors,  biased  or  unbiased.  Although 
the  quality  of  the  charges  was  excellent,  the  errors  are  surprisingly  large  for 
both  boostered  and  detonator  initiated  charges.  For  both  sets  of  data,  we  have 
omitted  from  the  graph  the  value  of  Rc  for  £  =  25. mm.  This  was  done  because  we 
are  not  sure  that  the  detonation  velocity  has  reached  a  steady  state  at  this  dis¬ 
tance  —  even  for  the  boostered  charge.  If  the  velocity  is  not  constant,  our 
method  of  analysis  breaks  down  --  see  the  derivations  of  Eg.  1  and  2.* 

Discarding  the  data  for  /  *  25.4  mm  leaves  only  2  points  to  establish  a  line 
for  radial  expansion  of  the  detonator  initiated  charges.  Tills  extrapolates  to  an 
effective  initiator  length  of  33  ran  for  RDX  (p0  *  1. 3  g/cc),  That  is,  the  front 
appeared  to  have  started  33  ram  behind  the  charge.  Certainly  more  data  on  curvature 
and  on  detonation  velocity  are  needed  in  this  case.  We  can  hazard  a  guess  that 
the  front  is  expanding  radially  out  to  // d  of  at  least  k.  5. 

The  results  for  pentolite  boostered  charges  of  RDX  ( Pq  =  1.3  g/cc)  are  also 
shewn  In  Fig.  9*  These  were  obtained  by  using  Eq.  1  (assuming  radial  expansion); 

errors  being  estimated  as  outlined  above.  A  straight  line,( - -)  fits  these 

results  fairly  well,  but  it  gives  an  effective  initiator  length  of  63  mm.  Sub¬ 
tracting  the  33  mm  obtained  above  from  the  63  mm  effective  initiator  length,  we 
have  30  mm  length  of  RDX  (1.3  g/cc)  equivalent  to  25,4  mm  pentolite  (3.56  g/cc). 
Inasmuch  as  the  ratio  of  the  impedances  (pentolite/RDX)  is  1.21,  this  seems  a 
reasonable  equivalence.  Apparently  we  do  not  have  enough  information  on  the 

*  The  two  points  at  £=  25.4  mm  [(//d)  =  O.5J  fall  well  below  the  curves  of 
Fig.  9.  The  measured  front  velocity  at  this  location  (see  Table  8)  supports 
the  suggestion  that  the  steady  state  velocity  had  not  yet  been  achieved. 

It  is  also  relevant  to  note  that  the  initiator  effect  is  still  dominant  at 
this  low  value  of  C^/d)  for  AF.  (See  Fig.  j). 
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boostered  system  or  on  the  detonator  initiated  system.  We  have,  however,  enough 
information  on  both  systems  so  that  we  can  design  electromagnetic  experiments  for 
measuring  particle  velocity  with  some  assurance  as  to  the  amount  of  curvature  to 
expect  in  the  detonation  front.  That  is,  the  front  in  RDX  (1.3  g/cc)  expands 
radially  out  to  at  least  (H/&)  —  .4. 5.  This  does  not  agree  with  the  behavior 
reported  for  1.2  g/cc  HDX  in  Fig  5.7a  of  Ref.  2.  The  difference  might  arise  from 
the  density  difference  (1.2  g/cc  vs  1.3  g/cc)j  from  a  difference  in  particle  size; 
or  from  a  difference  in  data  analysis. 

Tetryl 

The  curvature  of  the  detonation  front  in  tetryl  charges  has  been  reported 
previously*-.  A  35  ran  streak  camera  was  used  in  that  work  so  that  resolution  of 
both  time  and  distance  was  not  as  good  as  with  the  70  mm  camera  now  used.  One 
consequence  of  this  is  tnat  the  data  near  the  center  of  the  charge  were  of  low 
significance.  Thus  the  fit  to  Eq.  1  was  governed  mostly  by  the  data  over  the  last 
1.2  cm  on  each  edge  of  the  charge,  and  the  requirement  that  the  curve  pass  through 
the  points  t  a  0,  x  «  0.  These  data  have  been  re-examined;  they  were  fitted  to 
both  Eqs.  1  and  2.  The  main  change  is  that  the  results  for  the  200  mm  long 
charges  now  lie  close  to  the  radial  expansion  line  rather  than  above  it,  see 
Fig.  10.  This  is  probably  due  to  an  improvement  in  the  computer  cede.  Plotting 
the  experimental  data  and  the  values  computed  from  the  new  fit  clearly  demonstrate 
that  the  new  fit  is  superior  for  a  200  irm  long  charge.  Tetryl  ( p0  =  1.51  g/cc) 
exhibits  radial  expansion  out  to  200  ram^^/d)  =  h,  as  stated  previously. 

Summary  R^  vs  ^  Fig.  11  summarizes  the  R,,  vs  H  measurements  we  have  made  and 
reported.  Though  by  no  means  a  complete  survey,  it  strongly  indicates  that  any 
low  density  ( pQ  —1.0  g/cc)  explosive  will  rapidly  attain  a  constant  Rj,  at 
(£ /d)  >  1.5  whether  the  explosive  is  ideal  or  non-ideal.  High  density  conventional 
explosives,  e.g.,  tetryl  at  87$  TJ®,  exhibit  a  geometrical  spherical  expansion  of 
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the  detonation  front  at  Hfd  <  4. 0.  There  may  be  cases  intermediate  between  these 
two  extremes,  NQ  at  92$  TMD.  In  addition  to  oQ  ,  Rc  depends  on  the  nature  of 


the  material  and  on  its  particle  size. 


Estimates  of  Reaction  Zone  Length 


The  Wood  and  Kirkwood  curved  front  theory1*'  derives  the  relationship 

z  =  (s/3.5)(l-D/Di)  (4) 

where  z  is  the  reaction  zone  length,  s  is  the  radius  of  curvature  Rc,  and  the 
numeric  factor  is  somewhat  arbitrary.  For  H.E.  showing  constant  s  at  (i/d)  >  1. 5, z 
can  be  computed  from  our  present  data  although  the  results  for  NQ  and  TNT  will  be 
only  approximate.  For  expanding  fronts,  if  we  assume  z  constant,  then  s  and  D 
must  change  together.  Computation  of  z  then  requires  at  least  a  pair  of  (s,D) 
values ;  many  pairs  would  be  preferable.  Good  derived  z  values  for  tetryl,  and 
f'x  RDX  cannot  be  obtained  without  much  more  precise  measurement  of  both  D  and 
Rc  vs  but  an  estimate  of  the  size  of  z  can  be  made. 

Table  9  contains  the  results,  and  the  z  values  for  the  two  APs  are  plotted 
in  Fig.  12a;  Fig.  12b  shows  the  "a"  values  computed  from  the  original  curved  front 

c  C. 

theory'  and  reported  for  two  different  particle  size  APs  in  previous  work  . 

FLg.  12a  shows  first  that  the  density  has  only  a  small  absolute  effect  on  the  z 
value.  For  the  25  u.  AP,  vhatever  that  effect  is,  it  is  submerged  in  the  general 
scatter  of  the  results.  In  the  case  of  the  9  AP,  the  trend  is  quite  smooth: 
z  decreases  with  increasing  ^  to  a  minimum  at  1.2  g/cc  after  which  it  increases 
again.  The  minimum  occurs  at  about  the  density  where  the  D  vs  pQ  curve  departs 
from  linearity  (see  Fig.  4).  The  z  value  for  the  finer  material  is  less  than  that 
for  the  coarser,  but  the  ratio  value  is  0.24  to  0,19  whereas  the  6  ratio  of  9/25 
is  O.36.  Hence  the  ratio  of  these  nominal  reaction  zone  lengths,  z,  does  not 
equal  the  ratio  of  the  average  particle  sizes  as  did  the  ratio  of  the  nominal 
zone  lengths  a 
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In  Pig.  12b,  from  a  previous  report®,  both  the  25  H  and  the  10  11  material 

show  a  vs  pQ  curves  at  pQ  <  1. 0  g/cc  very  similar  to  the  z  vs  p^'  curve  of  t..e 

fine  AP  in  Pig.  12a.  Since  the  earlier  work  was  done,  we  have  found1^  that  for 

AP  at  pQ  >  1  g/cc  our  diameter  series,  terminating  at  a  =  7*  6  cm,  is  inadequate 

for  extrapolation  to  D*  values;  in  particular,  the  slope  of  the  D  vs  d'1  curve 

is  too  steep.  (Also  see  Appendix)  Because  "a"  is  directly  proportional  to  that 

slope,  the  a  values  at  p0  >1.0  g/cc  in  Pig.  12b  are  all  too  high,  and  because 

the  error  increases  with  the  density,  it  is  probably  responsible  for  the  sharp 

increase  in  *a”  indicated  at  higher  pQ  in  Pig.  12b.  In  the  region  where  this 

error  does  not  occur,  i.e. ,  pQ  <1.0  g/cc,  the  a/z  ratio  for  the  25  P1  AP  (N126 

in  both  cases )  is  about  3. 0.  The  9  H  AP  value  of  z  can  be  approximately  adjusted 

to  that  for  10  p,  by  assuming  proportionality  between  z  and  $ .  Then  the  a/ z 

ratio  for  10  ^  AP  is  3*9«  Pbr  AP  at  icw  f^,  the  a/z  ratio  seems  to  be  appreciably 

17 

less  than  the  value  of  7  reported  for  conventional  H.E.  1 . 

The  value  of  z  for  the  9  H  AP  at  1.3  g/cc  (1.10  mm)  is  quite  close  to  a  pre¬ 
liminary  measurement  of  the  zone  by  the  electro-magnetic  method  (1.1  to  1.7  nun  at 
1.29  g/cc).  The  z  value  for  low  p0  TNT  compares  fairly  well  with  measurements1® 
(1.2  vs  1.6  mm)  but  that  for  moderate  pQ  KDX  does  not  (0.2  mm  vs  0.6  to  0.9  mm). 
In  both  these  cases,  however,  the  agreement  or  disagreement  could  result  from 
very  small  errors  in  the  value  of  D  used.  For  example,  Table  8  lists  three  accep¬ 
table  experimental  values  of  D  for  EDX;  each  is  within  0. 6$  of  D^.  For  spherical 
expansion  of  the  detonation  front,  Rq  and  D  are  both  assumed  to  increase  with 
(//d)  until  ideal  values  (  “and  D^)  are  reached.  Hence  the  two  values  of  D  < 
of  Table  8  were  chosen  to  compute  z  (Table  9);  the  apparent  agreement  for  the  two 
results  is  probably  fortuitous  and  depends  on  the  differer  -,f  0. 02  mm/nsec  in 
D  which  is  well  below  our  experimental  resolution.  Future  work  will  include  addi¬ 
tional  investigation  of  experimental  measurement  and  methods  of  computation  of 
reaction  zone  lengths. 
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SUMMARY  AND  CONCLUSIONS 

The  present  study  of  curvature  of  the  detonation  front  was  carried  out  at 
constant  diameter  ( ~  50  mm).  Hence  no  information  on  the  effect  of  d  on  Rg  was 
obtained  although  d  is  considered  one  of  the  variables  which  determine  Rc. 
Information  obtained  in  the  study  follows: 

1.  Non-planar  detonation  fronts  are  well  approximated  by  a  spherical 
shape;  edge  effects  are  small. 

2.  Rc  is  a  function  of  the  chemical  material,  its  particle  size  3,  and 
its  Rg  increases  with  decreasing  5.  Rc  may  increase  with  in¬ 
creasing  P  (e.g. ,  NQ)  or  decrease  with  increasing  p  (e.g.,  AP). 

0  0 

The  effect  of  pQ  nay  be  associated  with  Group  1  or  Group  2  classi¬ 
fication. 

3.  Any  charge  having  pQ  <  1.0  g/cc  appears  to  have  a  constant  R^ 
rapidly  attained,  and  independent  of  the  type  of  initiation.  For 
point  initiation,  constant  Rc  is  established  at  (£/&)  >  1.5* 

4.  A P  exhibits  rapidly  attained  constant  curvature  of  the  detonation 
front  over  the  density  range  in  which  it  Is  detonable,  whereas 
tetryl  (and  presumably  other  conventional  H.35. )  at  high  p0  (>87#  ©ffl) 
exhibits  a  geometrical  spherical  expansion  with  i/d.. 

5.  NQ  may  be  excluded  from  "conventional"  H.E.  since  at  92$  3MD 
it  appears  to  exhibit  a  behavior  intermediate  between  that  of 
AP  and  that  of  tetryl  at  87$  TMD. 

6.  F!rom  the  standpoint  of  failure  theory,  the  behavior  of  subcritical 
charges  is  of  much  interest.  Strongly  boostered,  subcritical 
charges  show  failure  starting  with  quenching  of  the  reaction  at 
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t'ie  charge  surface  and  working  toward  the  axis.  There  is  no 
marked  increase  in  Rc  during  this  process,  nor  does  Rc  seem 
to  decrease  to  0.5  d  at  d  *  dc. 

7.  The  method  of  correcting  our  measured  D  values^  is  not  applicable 
to  any  charge  at  r<0  ~1.0  g/cc  or  to  any  charge  of  AP.  It  does 
seem  to  be  applicable  to  most  high  density  organic  H.E.  although 
it  is  only  partially  so  to  NQ  at  92$  M3.  Consequently  we  can 
use  this  correction  only  in  the  special  cases  where  we  know  or 
have  shown  by  test  that  it  is  valid. 

8.  In  the  electromagnetic  measurement  of  particle  velocity  a  plane 
wave  booster  cannot  produce  a  steady  state  planar  detonation  in 
AP  charges  or  (probably)  any  other  charges  of  p0  ~1.0  g/cc.  For 
these  charges,  the  constant  D  and  Rc  will  be  established  at 

(//d)  >1*5.  (We  have  previously  assumed  that  (//d)  >  3  is  necessary). 
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FIG.  3  TWO  MODELS  FOR  EXPANSION  OF  DETONATION  FRONT  IN  CYLINDRICAL  CHARGE. 
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FIG.  5  FAILING 
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FIG.  6  Rc  VS  p  CURVES  FOR  9fi  AND  25fi  AP  AT  d~50mm. 
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FIG.  7  RADIUS  OF  CURVATURE  VS  LENGTH  OF  AP  CHARGES  OF 
N-138  AND  N-126,  d  =  50.8mm 
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FIG.  10  Rc  VS  LENGTH  OF  TETRYL  CHARGES,  d  =  50.8  mm 
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(a)  NOMINAL  REACTION  ZONE  LENGTHS  ,  z,  FOR  AP  AS  A  FUNCTION  OF  DENSITY  &  PARTICLE  SIZE 


LOADING  DENSITY  PQ  (g/cc) 

(b)  NOMINAL  REACTION  ZONE  LENGTHS,  FOR  AP  AS  A  FUNCTION  OF  DENSITY  &  PARTICLE  SIZE 

FIG.  12  NOMINAL  REACTION  ZONE  LENGTHS  FOR  AP  AS  A  FUNCTION  OF 
DENSITY  AND  PARTICLE  SIZE. 
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Table  1 

DETONATION  VELOCITY  DATA  FOR  AP»S 
P0  d  D 

.Srot  Nq^ _ (g/sg-1  _ (afflX.  (nwjxsec  )  = _ Comment - 


AP(N138)  8 

49t* 

651 

0.627h 

50.8 

2.452 

654 

0.633h 

50.8 

2.502 

704-1 

0.920H 

50. 6a 

3.286 

705-1 

0.920H 

50. 6a 

3.302 

706-1 

0.970H 

50.6s 

3.459 

710-1 

0.970H 

50.6s 

3.394 

652 

1.02CH 

50.3 

3.558 

655 

1.020H 

50.3 

3.590 

697 

1.1901 

49.9 

3.985 

698 

1.2991 

49.9 

4.008 

708 

1.3451 

49.9 

4.190 

752 

1.354i 

49.9 

4.340 

709 

1. 384i 

49.8 

4.317 

699 

1.4111 

49.9 

F 

700 

1.4871 

49.9 

F 

AP(N136)b 

8.4n 

602 

1.074H 

50.8 

3.713 

590 

1.2901 

50.8 

4.188 

AP(N126)c  25  tt 


693 

0. 901h 

2.134 

Poor  quality  charge  and  trace. 

703 

l.OOlh 

2.393 

694 

1.008h 

2.401 

714 

l.lOlh 

2.586 

695 

l.lOlh 

2.683 

689 

1.206i 

2.794 

691 

1.246i 

2.859 

690 

1.3021 

2.8$0 

702 

1.346i 

2.612 

Marginal  for  detonation. 

692 

1.3861 

(2.475 )F 

Charge  length  196.9  mmj 

all  others  228.6  mm 

b.  See  Appendix  for  additional  data 

c.  d  =  50.8  min,  all  charges 
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Table  2 

SUBCRmCAL  REACTIONS  IN  9H  AP(N138) 


Ax/At  (nim/(Asec) 


Mid-Pin  Location 
nan 

P0  =  1-59 
g/cc 

P0  =  1.635 

_ sZsi l _ 

25.4 

4.47 

4.58 

50.8 

4.39 

4.62 

76.2 

4.43 

4.4o 

101.6 

4.37 

4.15 

127.0 

4.17 

3.50 

152.4 

3.64 

2.70 

177.8 

3.42 

1.83 

203.2 

2.43 

1.37 

228.6 

1.93 

No  reading.  Last 
two  inches  of  tube 
recovered  undamaged 
except  split  in  two, 

Charges  in  seamless 

steel  tubes  3*65  cm  I.D. , 

4.76  cm  O.D. , 

30.5  cm  length. 

Booster  of  50/5O  pentolite  (1.56  g/cc),  5.08 

cm  diam  x  5.08 

cm  long 

Pins  spaced  at  25.4 

mm  intervals  starting  12.7  mm  from 

booster  surface 
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Table  3 

RADIUS  OF  CURVATURE  AS  A  FUNCTION 
OF  DENSITY  IN  AP 


R„  (±  $  error),  mm 


Shot  No.  Po  (g/cc)  Whole  Front  Middle  of  Front 

9  h  AP,(H38) 


725 

0.922H* 

70.3 

1.57) 

97.8 

[4.47) 

734 

1.020H* 

7  6.6 

1.23] 

86.4 

2.01) 

666 

1.075H 

75.8 

1.52 

81.6 

1.69) 

667 

1.075H 

74.5 

1.57! 

88.8 

*2.34) 

668 

1.075H 

70.91 

1.50 

88.71 

;$.78 

713 

1.1911* 

66. 3 

!l.26) 

75.1 

3.12} 

711 

1.298i* 

67.0 

1.35 

76.5 

1.32) 

748 

1.3471* 

63.2 

1.22] 

71.4 

1.75) 

735 

1.376i* 

64.9 

1.14] 

76.3 

[1.52^ 

723 

724 
656 
662 
6t4 
677 
676 
685 


25  n  AP(N126) 


l.OOlh 

51.01 

1.17 

56.81 

(i.72) 

l.lOlh 

55.21 

0.85) 

62.31 

1.38) 

1.191ia 

60. 4l 

1.60) 

70.11 

1*37 

1.258i 

53.1(0.96) 

63.ll 

ta.4l) 

1.273ib 

52.01 

(0.64) 

54.91 

(1.71) 

1.3381 

43.51 

(0.86) 

50.31 

(1.07) 

1.3391 

44.71 

■0.82) 

49.4(1.68) 

1. 3'flic 

46.81 

(1.01) 

54.7(0.69) 

*  d  =  49.9  urn.  For  all  other  changes,  d  =>  50.8  mm 

a.  203.2  mm  long 

b.  457*2  mm  long.  All  others  228.6  mm  long 

c.  This  density  may  be  above  pcj  see  Fig.  4. 
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Table  4 


RADIUS  OF  CURVATURE 
va 

CHARGE  LENGTH  FOR  AP 


Rg  (±  j,  error),  mm 


Shot  No. 

PQ  (g/cc> 

0(mm) 

Apparent  D 
nun/p.sec 

From  Eq.  2 

From  Eq. 

742 

1.191 

9  H  (N138) 

25.4  4.27 

74.2(1.12) 

74.5(2.2) 

743 

1.191 

25.4 

4.27 

72.2(1.12) 

72.1(2.5) 

726 

1.193 

76.2 

4.08 

67.8(0.95) 

68.5(1.7) 

727 

1.191 

101.2 

3.99 

69.7(1.25) 

68.5(1.4) 

728 

1.200 

127.0 

4.05a 

67.7(1.26) 

68.3(3.8) 

729 

1.191 

177.8 

3.99 

68.8(1.19) 

- 

713 

1.191 

228.6 

3.99 

66.3(1.26) 

- 

744 

1.342 

25.4 

4.36 

72.6(1.29) 

73.1(2.0) 

745 

1.345 

7o.2 

4.36 

64.0(1.14) 

65.3(2.4) 

746 

1.358 

127.0 

4.46 

63.6(l»4i) 

64.4(2.3) 

747 

1.348 

177-8 

4.36 

63.5(1.10) 

63.2(1.22) 

- 

748 

1.347 

228.6 

4.36 

- 

75  6 

1.296 

25  u  (N126) 

25.4  3.53 

71.9(1.71) 

757 

I.296 

76.2 

3.62 

49.0(0.78) 

- 

b 

I.296 

228.6 

- 

49.3  - 

- 

1 


*  All  charges  prepared  in  isostatic  press;  d  =  49.9,  50.8  mm  for  IO.38  and 
N126,  respectively. 

a.  Correct  to  4.02  for  density  of  1.191 

b.  Read  from  Fig.  6. 
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Table  5 

FRONT  CURVATURE  FOUND  IN  SUBDETONATION 
REACTION  OF  AP  CHARGES 


Shot  No. 

P„(g/cc) 

/(mm)  D(cm/^sec) 

R„  (±  $>  error)  mm 

Comment 

753 

754 

1.4lli 

1.409i 

9  |i  AP(N138) 

76.2  (4.385) 

25.4  4.47* 

70.7(1.17) 

Obvious  curvature 
Trace  covers 

755 

1.409i 

76.2  4.38* 

61.0(0.96) 

charge  diam 

Trace  covers 

683 

1.399i 

26  \i  AP(N126) 

152.4 

Not  read 

charge  diam 

Trace  covers 

682 

1.4o8i 

177.8 

Not  read 

charge  diam 

Very  faint  trace 

730 

1.4351 

76.2  2.44* 

59.7(0.81) 

covers  only 
central  third 
of  diam 

Trace  covers 

731 

1.435i 

114.3  2.44* 

48.9(1.05) 

charge  diam 
Trace  covers 

_ _ 1  _  . 

only  9C$  diam} 
see  Fig.  2a 


* 


Shock  velocity  used  in  computing  Rc 


RADIUS  OF  CURVATURE  vs  CHARGE  LENGTH 
FOR  NQ-i  (X5^7) 
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Table  7 

RADIUS  OF  CURVATURE  MEASUREMENTS  ON  TNT,  X517a 
AT  p0  =  l.C6l  g/cc 


782  76.2  95. 6C 

783  76.2  86.6  Av.  92.04  {1/ d)  =  1.82 

784  127*0  88.1 

785  228.6  97.3 

All  charges  were  5G8  mm  diam  and  were  prepared  in  the 
hydraulic  press 

a.  Rotap  sieve  analysis  on  lOOg  sample: 

Retained  on 


Screen  No. 

wt,  .(g) 

30(590  n) 

0 

50(297  n) 

5.20 

100(149  \i) 

67.20 

140(105  ;i) 

26.50 

Pan 

0.10 

99.00 

Average  particle  size  ca.  200  |J. 

b.  D  of  4.8  mm/p.sec  (about  8$  below  ideal)  used  to  convert 
measured  time  to  distance.  The  larger  correction  to  Dj^ 
than  that  of  Table  6  for  NQ-/  was  an  estimate  for  the 
different  charge  preparation. 

c.  Bottom  part  of  this  trace  was  blocked  out;  cause  unknown. 
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Table  9 


NOMINAL  REACTION  ZONE  LENGTHS  (z) 
COMPUTED  FROM  EQ.  4 


Velocity  (am/ usee) 


D 

(d  =  50.8) 

Fig.  4 

D-i** 

Curvature 

s(xm) 

FigsTb  &  11 

Zone 

z(mui) 
Eq.  4 

1.45 

0.9 

3.23 

9  U  AP 

3.464 

75.2 

1.0 

3.49 

3.722 

72.7 

1.30 

X#  X 

3.76 

3.980 

70.3 

1.11 

1.2 

4.01 

4.237 

67.8 

1.04 

1.3 

4.23 

4.495 

65.4 

1.10 

1. 35 

4.30 

4.624 

64.1 

1.28 

1.0 

2.38 

25  u  AP 

3.722 

50.8 

5.23 

1.1 

2.62 

3.980 

55.0 

5.36 

1.2 

2.81 

4.237 

59.3 

5.69 

1.  3 

2.83 

4.495 

48.8 

5.15 

1.35 

2.68 

4.624 

43.5 

5.21 

1.001 

5.2* 

m-i 

5.455 

94 

1.3 

1.08l 

5.07* 

TNT 

5.318 

92 

1.2 

1.30 

7.01* 

RDX 

7.05 

130.5 

0.21 

1.30 

7.03* 

7.05 

270.9 

0.22 

*Better  values  needed  --  especially  for  TNT  for  which  (Di  -  D)  was  equated 
to  ttu  'j  found  for  N Q-/  (l  g/cc). 

^Respective  references  are:  A ?r  (12);  NQ,  (7);  TNT  and  RDX,  (l6). 
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APPENDIX 

Additional  Data  on  AP*  s 

When  NCS*  began  preparing  batches  of  AP  which  were  obviously  finer  than  the 
10  (i  average  particle  size  of  our  initial  work®,  one  more  attempt  was  made  to 
obtain  higher  density  D  vs  d  data  (with  d  <  7.6  cm**)  which  extrapolate  to  the 
Dj_  vs  pQ  curve  derived  at  pQ  «  1,0  g/cc1^.  Table  A1  contains  the  D  vs  d  data 
for  8.4  H  AP  (H136)  at  1.074  and  1.30  g/cc.  Fig.  A1  shows  the  extrapolation 
and  Table  A1  compares  the  values  so  obtained  with  Dj  given  by  the  analytical 
expression  of  Ref.  (12).  When  the  series  is  terminated  at  d  *  7.6  cm,  the  extra¬ 
polated  value  is  greater  than  the  analytical  by  4$  at  p0  *  1.074  g/cc  and  by  12 ‘jo 
at  po  =  1.30  g/cc.  These  are  about  the  size  differences  found  in  the  previous 
work.  It  should  be  noted  that  some  of  the  difference  found  can  be  attributed 
to  omitting  the  correction  (now  shown  inapplicable  to  AP  charges)  to  the  measured 
D  values.  If  such  a  correction  were  used,  the  difference  (at  p0  *  1.074  g/cc) 
would  have  been  reduced  to  2f>  but  would  not  have  been  eliminated. 

An  apparent  aging  effect  (decrease  of  D)  in  AP  was  reported  earlier12. 

This  has  been  confirmed  by  the  results  on  25  U  AP  reported  in  the  text  as  well 
as  those  shown  in  Table  A2  for  an  8.8  (i  AP.  The  difference  found  is  small  but 
consistent,  a  decrease  of  about  5$  in  D  with  age.  Differences  with  age  can  some¬ 
times  be  seen  in  the  AP.  Thus  the  N139  of  Table  A2,  although  it  had  been  opened 
only  in  a  dry  box  or  heated  in  a  vacuum  oven,  showed  some  caking  with  age  and  had 
to  be  rolled  before  making  charges  721  and  722.  Because  we  cannot  handle  this 
fine  AP  without  causing  some  agglomeration,  henceforth  charges  will  be  made  from 
it  without  any  oven  treatment  of  the  AP. 

*Naval  Ordnance  Station,  Indian  Head,  Maryland. 

**Charges  small  enough  to  do  no  damage  to  our  bombproof  equipment. 
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Table  A1 

D  VS  d  DATA  FOR  AP(KL36),  8.4ia 


Shot  No. 


599 

600 
602 
603 


587 

588 

590 

591 


d  ('em  ) 

'(g/fcc) 

D(mm/hsec) 

2.54 

1.072H 

3.026 

3.495 

1.073H 

3.525 

5.08 

1.074k 

3.713 

7.62 

1.073H 

3.837 

2.54 

1. 304i 

(2.749)F 

3.495 

1.3021 

3.787* 

5.08 

1.2901 

4.188 

7.62 

1.2951 

4.461 

D-i  tnm/u«oc 
Extrap.  Ref  (12) 


4.06  3.91 


5.04  4.49 


*  Poor  record. 
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Table  A2 

AGING  EFFECT  ON  AP(N139),  8,8  \i 


Shot  No. 

d(mm) 

A^(s/cc) 

Date 

Dfmm/nsec^ 

672 

673 

50.6 

50.6 

I-* 

II 

Pm 

0  0 

3.466 

3.508 

Av.  3.49 

721 

722 

50.8 

50.8 

1.001H 

1.001H 

10/25/68 

10/25/68 

3.315 

3.292 

Av.  3.30 


Charge  Length:  228.6  mm 
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